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Abstract 

In this paper a quantum mechanical phase space picture is constructed for coarse- 
grained free quantum fields in an inflationary Universe. The appropriate stochastic 
quantum Liouville equation is derived. Explicit solutions for the phase space quantum 
distribution function are found for the cases of power law and exponential expansions. 
The expectation values of dynamical variables with respect to these solutions are com- 
pared to the corresponding cutoff regularized field theoretic results (we do not restrict 
ourselves only to ($^)). Fair agreement is found provided the coarse- graining scale is 
kept within certain limits. By focusing on the full phase space distribution function 
rather than a reduced distribution it is shown that the thermodynamic interpretation 
of the stochastic formalism faces several difficulties (e.g., there is no fluctuation- 
dissipation theorem). The coarse-graining docs not guarantee an automatic classical 
limit as quantum correlations turn out to be crucial in order to get results consistent 
with standard quantum field theory. Therefore, the method does not by itself consti- 
tute an explanation of the quantum to classical transition in the early Universe. In 
particular, we argue that the stochastic equations do not lead to decoherence. 
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I. Introduction 

The paradigm of stocliastic inflation, flrst introduced explicitly by Starobinsky 0, 
has recently become popular as a means of investigating various features of inflation. 
Some studies using this formalism are those of density perturbations from inflation 
0, the very large scale structure of the Universe 0, "eternal inflation" power 
law inflation 00, and speculations regarding the relationship of this formalism to 
quantum cosmology ||^ (this list is by no means exhaustive). 

It must be admitted, however, that there is still no iron-clad justification for the 
systematics of the method nor, for that matter, a clear-cut interpretational scheme. 
The claim at issue is that the infrared behavior of massless or small mass quan- 
tized scalar fields in an infiationary Universe can be described in terms of a real 
time classical random process. The source of the noise is taken to be large scale 
quantum fiuctuations which are continuously generated in an infiationary Universe 
by red-shifting of the ultraviolet sector. A key question here is: can these quantum 
fiuctuations be treated as being classical? 

These fundamental issues have been considered previously for free fields ||^ how- 
ever the situation for interacting fields is not clear, and it is not obvious how far, if 
it all, any of the present "derivations" are correct In this paper we leave aside 
for the moment the problem of interacting fields and attempt a further clarification 
of the issues addressed in Ref. p. To do so we will derive a phase space quantum 
master equation for a quantum distribution function (the Wigner function) and study 
its solutions. In our stochastic approach averages with respect to this distribution 
function are supposed to reproduce quantum field theoretic expectation values. A 
study of the solution itself is supposed to enable one to judge the "classicahty" of 
each physical situation. The new method is distinct from the conventional approach 
(where one takes as given a classical Langevin equation), and enables the inclusion 
of crucial quantum correlations that have been unjustifiably neglected in the past. 

A subtle and important aspect of the Wigner distribution function is the fact that 
essential quantum features are hidden in quantum correlation "cross-terms" that dis- 
appear when one integrates over any one of the phase space variables to produce 
a one-variable (necessarily positive definite) distribution function. Such a reduced 
distribution is essentially useless as a diagnostic tool for studying quantum correla- 
tions in phase space (as will be seen forcefully in this paper). Unfortunately it is on 
precisely such objects that attention has been focused till now. Here, with the full 
Wigner function at hand we will be able to go much further with regard to clarifying 
the physics behind the stochastic approach. It has been noticed previously that 
the reduced distribution for a massive scalar field in de Sitter space has, at late times. 
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an intriguing thermodynamic interpretation: it corresponds to a Boltzmann distri- 
bution at the Gibbons-Hawking temperature. However, the full distribution found in 
this paper does not have a thermal form even though the reduced distribution is the 
same. This can be traced directly to the fact that quantum correlations have not been 
neglected; indeed they are every bit as important as the remaining contributions. We 
will go more deeply into this question in Sec. V. 

It has long been appreciated that the stochastic approach probes the infrared 
sector of the relevant field theory. The length scale is set by a certain parameter e 
which is usually taken to be small (i.e., attention is confined to length scales much 
larger than the time dependent horizon length). Assuming inflation began at a finite 
time in the past, one cannot take e to be arbitrarily small independent of the time scale 
of interest: a small e is consistent only with "late times." It is also well known that 
the quantum theory of free fields in an inflationary spacetime has a nontrivial infrared 
sector and that a simple (though certainly not rigorous) way to calculate expectation 
values of fleld variables is to set an upper momentum cutoff at the Hubble scale 
(corresponding to e ~ 1). The stochastic picture conflicts with these fleld theoretic 
results unless e is small; this is due to the fact that in the stochastic approach one 
focuses essentially on the zero mode and attempts to include inhomogeneities only 
through a noise term. Though this approximation appears to be quite drastic, we will 
show that the stochastic calculations even for e ~ 1 are never too far from the naive 
field theoretic results. (Unlike some previous work, our formalism does not restrict 
the value of e.) 

Previous work in stochastic infiation has concentrated mainly on the quantity 
($^). In this paper we extend the method to compute {^n + n^) and (tt^). Such 
quantum averages would be needed if one wished to compute the expectation value 
of the stress tensor. (While for a small mass field, and an exponential infiation, the 
dominant contribution to to (T^i/) is from terms oc ($^), it is important to check if 
the calculation of the other terms is trustworthy.) Earlier approaches to stochastic 
infiation implement approximations which led to incorrect values for these quantities. 
Indeed precisely these approximations formed the basis of some arguments about 
stochastic infiation leading to an automatic classical limit for the quantum field theory. 
We will argue against any such result in Sec. V. 

An interesting (and somewhat uncomfortable) feature of the quantum phase space 
distribution found in this paper is that, in some cases, it depends quite strongly on 
e and is indeed singular in the limit e — 0. Therefore, while it is true that ($^) (as 
well as the reduced distribution for the fleld alone) may be independent of e for small 
values of e, this parameter does not drop out of the physics. A finite value of e is 



4 



necessary for the distribution function to exist; this is true even for a massive field 
in an exponentially expanding Universe where to leading order all quadratic phase 
space expectation values are independent of e. 

An important issue that seems to have received insufficient attention in the stochas- 
tic inflation literature is a discussion of the role of initial conditions. Massless theories 
in inflationary spacetimes suffer from infrared divergences. Typically these diver- 
gences are "fixed" by assuming that infiation began a finite time in the past and 
thereby modifying the infrared structure of the quantum state of the field. Expec- 
tation values then have two contributions: one each from the pre-infiationary and 
infiationary sectors. One can show that the pre-infiationary contribution falls rapidly 
with time and can always be neglected compared to the infiationary one (see the Ap- 
pendix). In the stochastic paradigm there are also two contributions to expectation 
values: a systematic piece arising from the dynamical evolution of an initial condition 
and a stochastic piece due entirely to the noise source. We will show that while at 
late times, and for arbitrary initial conditions, the second piece always dominates the 
first, this is not true at early times. While the matching of the quantum states in the 
pre-infiationary and infiationary regimes at the onset of infiation can also provide an 
initial condition for the stochastic method, the time dependences of the systematic 
contribution in the stochastic method do not always match the time dependences of 
the pre-infiationary contribution in the field theoretic calculation. This fact coupled 
with the small e restriction might limit the application of stochastic techniques in 
accurately studying the onset of infiation. The use of the method for studying phase 
transitions in the early Universe should also be approached with some caution [|10| . 
(Of course all this is not a serious problem if one is only interested in late time results.) 

The attempt in this paper is to push the formalism of stochastic infiation as hard 
as possible in simple examples: we find that some of the appealing original results 
no longer appear as compelling as at first sight. However, it is still a remarkable fact 
that a simple stochastic model suffices to (almost) correctly calculate field theoretic 
expectation values and further that the essentially nonstationary phase space distri- 
bution nevertheless yields a thermal (or "random walk" ) distribution for the reduced 
distribution function. Whether this has a deep significance is unfortunately not clear. 

The organization of the paper is as follows: In Sec. II we derive the appropriate 
stochastic quantum Liouville equation for the coarse-grained field using the phase 
space formulation of quantum mechanics and obtain the general solution. In Sec. Ill 
we apply these results to the case of an infiationary expansion; power law expansions 
are dealt with in Sec. IV. The existence of classical stochastic interpretations is 
discussed in Sec. V via a study of the solutions of the quantum Liouville equation. 
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We conclude with Sec. VI where the results are reviewed and future directions for 
research are suggested. The quantum field theoretic derivations of the results obtained 
via the stochastic approach are given in an Appendix. 
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II. The Stochastic Quantum Liouville Equation 

In this section we will set up a formalism to study the evolution of coarse-grained 
free scalar fields in a spatially flat inflationary Friedmann- Robertson- Walker Universe. 
We will work under the "test field" assumption, i.e., the contribution to the stress 
tensor from the field is taken to be small compared to that of the matter driving the 
expansion. All of our results will therefore not be applicable to an infiaton field but 
some may indeed be extended to that case. 

The line element for the spacetime is 

ds^ = —dt^ + a{tf'dx-dx (1) 
= S{riY (-drf + dx ■ dx) , (2) 

where in terms of the cosmic time t, the conformal time 

A massive minimally coupled scalar field has the Lagrangian 

£($, $,^) = -^v^ (9"'^,,^,u + m2$2) , (4) 
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which, with the metric choice (|^), reduces to 

The overdot represents differentiation with respect to the conformal time. In terms 
of the "conformal field," defined via the time dependent canonical transformation, 

X = S^, (6) 

and modulo an integration by parts, the Lagrangian becomes 

(7) 



-x^ + x.iX.i + ( 5W - ^1 x^ 



s. 



One advantage of working with the conformal field and the Lagrangian (|^) is that the 
equation of motion for the field 

X - V'X + (^S^m^ - 1^ X = (8) 

does not have the first derivative in time "Hubble damping" term found in the equa- 
tion of motion for the original field 

<!> + 2|<i> - V^* + 5W<I> = 0, (9) 



and that the canonical momentum 

TT,, = 

dX 



dC . 

t:^= =X (10) 



is of the usual flat space form. 

The Hamiltonian corresponding to (|^ is 



H{X, '^x) = \j dx 



-K^ + X^iX^i + (sW - I) 



which is that for a free field in flat spacetime with a time dependent mass. The 

Hamiltonian equations of motion are 

5H , , 

X = ^ = vr,, (12) 



The form of the Hamiltonian (pUf ), though "canonical," is hardly unique; if we 
had worked with some other choice of time and field it would have been "natural" 
to consider a different description in terms of a different Hamiltonian. At the level 
of a classical treatment, and even at the level of quantum dynamics, this difference 
is largely irrelevant. However, there is an aspect of the quantum treatment where 
such a difference does indeed matter: The two Hamiltonian descriptions will, upon 
canonical quantization, lead to inequivalent descriptions in terms of different vacua. 
Furthermore, there are well known difficulties if one chooses to select the time depen- 
dent ground state of the Hamiltonian as the "instantaneous diagonalization" vacuum 



11 1 . These problems will be of no concern to us as in our case the choice of quantum 



state will be an independently defined adiabatic vacuum (details will be given later). 
(We note in passing that for spatially fiat FRW models, Weiss has shown |12[ that 
with the specific Hamiltonian (|l^) the instantaneous diagonalization approach can be 
made consistent with a "mode quantization" for certain special choices of the latter.) 

In our case there is a conceptually important consequence of (|lTj) being the chosen 
Hamiltonian. In the stochastic inflation literature there appears to be a tendency of 
interpreting the Hubble damping term in (|^) as being of a truly dissipative nature. 
This runs the risk of repeating an old error in quantum mechanics: the confusion of a 
time dependent mass with true damping Working with ( |Tll) and the associated 
equation of motion (^ manifestly eliminates the possibility of such misinterpretations. 

The scalar field is now quantized in the standard manner first we introduce 
the modes 



M^,V) = S{ri)^j;ix,ri) (14) 
;MV) (15) 



(27r)3/2 



where Xk{ri) is a solution of 

Xk + uJ^^Xk = (16) 
with ujk the oscillator "frequency," defined via 

u,^ I^S^m^ - I) . (17) 

The annihilation and creation operators with respect to these modes, which satisfy 
the commutation relations, 

i^k.ai] = 5{k-k'), (18) 

[H^H'] = (19) 
are then used to build the field operator 

X{x,r]) = S{7])^x,r]) (20) 



dk 



«fc^fc(^>^)+4^K^'^)l- (21) 



As is well known the annihilation and creation operators are not uniquely specified by 
(p!8|) and (|19]). Further restrictions are needed to fix these operators and thereby to 
uniquely specify the "vacuum" state annihilated by a^. We will turn to these issues 
shortly. 

A seemingly generic feature of inflationary spacetimes is the "destabilization" of 



massless scalar fields [|T5] ||T6| due in part to infrared divergences . To render the 
quantum state infrared finite one assumes a benign Robertson- Walker expansion in 
which there is no infrared divergent adiabatic vacuum prior to inflation (a radiation 
dominated Universe, for example). The quantum state in the inflationary phase is 
matched to an infrared finite quantum state (e.g., the conformal vacuum) at the time 
when inflation takes over from the previous epoch. It is then possible to show that 
the new state is always free from infrared divergences |T^. The key result, however, 
is that the expectation value (<l>^) in the infrared finite (and ultraviolet regulated) 
state starts to grow at the onset of infiation; for power law infiation ($^) grows to an 
asymptotic constant value, whereas in the case of an exponential expansion it grows 
linearly with cosmic time without any upper limit. These otherwise puzzling results 
have a natural interpretation within the framework of stochastic infiation p|. 

Crudely speaking, "destabilization" occurs when, with m = 0, the cu^^ term in 
(p!6|) goes negative with the passage of time, at ever higher values of k. Modes at long 
wavelengths behave as amplitudes for upside down harmonic oscillators (with time 
dependent "frequencies" ) and due to the infiationary expansion there is a continuous 
fiow of short wavelength modes into this unstable infrared sector. Therefore we focus 
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attention on the long wavelength modes, i.e., those with k"^ < S/S by defining the 
coarse-grained quantum field 

(22) 



Xiix, r]) = J dk 6{ks - k) d^X^{x, r]) + dix*:{x, r]) 
and the corresponding coarse-grained momentum 



TTiix, 7]) = I dk e{ks - k) d^X^{x, rj) + dtxHx, rj) 



(23) 



For the moment we will leave the upper cutoff ks unspecified beyond the fact that it is 
set by S/S (however, it is important to remember that this cutoff is time dependent). 
The corresponding short wavelength fields Xs and are defined by 



X 



Xr + X, 



5; 



(24) 
(25) 



The Heisenberg operators X and vr^ satisfy the classical Hamiltonian equations of 
motion. Substituting ( pi] ) and ( P5| ) in (^) and (|TB|) we find 



Xl 



VXl 



S 



c 

2 ) 



where 



(f, 7]) = ks j dk 5{k- ks) a,^X^{x, t]) + atxi(f , r]) 
F2{x, 1]) = ks dk 6{k- ks) [dj;Xj;{x, t]) + atx|(f , r]) 



(26) 
(27) 

(28) 
(29) 



The new terms and F2 arise simply because ks is time dependent. These terms 
represent the infiow of short wavelength modes into the infrared "condensate." (It is 
important to note that this contribution exists even for free fields.) 

In order to proceed further we have to decide which quantum state the field is in 
during the infiationary phase. It is known that the adiabatic vacuum suffers from in- 
frared divergences; to produce states free of such divergences one usually modifies the 
long wavelength mode structure of the quantum state (i.e., long compared to the time 
dependent horizon length at the onset of infiation) but leaves the short wavelength 
structure the same. This implies that the quantum state for computing expectation 
values of and F2y of various powers of these operators, is the adiabatic vacuum. 
Other choices are possible when describing different physical situations, for example, 
thermal states have been considered in Ref. and more general vacuum states in 
Ref. m. 
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Eventually we will deal specifically with an exponential expansion and with power 
law inflation (i.e., where the radius of the Universe goes as a power, greater than 
one, of the cosmic time). For such a Robertson- Walker Universe we will assume the 
quantum state for the short wavelength modes to be the adiabatic vacuum (which 



reduces to the Bunch-Davies vacuum [T^ for de Sitter space). For the moment, 
though, all that is relevant is that the chosen state be annihilated by the operator 
of (BTD, so that, 



It is also straightforward to compute that 



0, 
0. 



K[{xi,rii)F^{x2,ri2)) = 2Bij{xi,X2,r]i)S{rii - V2) 



where, with R = \xi — X2\, 

Bu{xi,X2,r]i) 
Bi2{xi,X2,r]i) 
B2i{xi,X2,r]i) 

B22iXi,X2,r]l) 



^k- 
^k- 



ks 
ks 
ks 



sin ksR 2 

sin ksR 

ksR 
sin ksR 

ksR 
sin ksR 

ksR 



(30) 
(31) 



(32) 

(33) 
(34) 
(35) 
(36) 



It is at this point that a stochastic interpretation suggests itself. The quantum ex- 
pectation value may be regarded as an averaging bracket for the white (albeit non- 
stationary) "noise" operators and Since we are dealing with a free theory it 
is trivial to verify that the higher moments of these operators are those appropriate 
for Gaussian noise. The fact that the noise is white stems from the theta function 
cutoff in momentum space. Other cutoffs are certainly acceptable, however, they will 
lead to the noises being colored and unnecessarily complicate the derivation of the 
phase space picture. We emphasize that physical results do not depend strongly on 
this choice. 

The "diffusion matrix" Bij has two curious features: it is complex (albeit Hermi- 
tian) and singular. We will show later that as far as the stochastic quantum Liouville 
equation is concerned what is really relevant is Bij + Bji which not only is necessarily 
real but also has a nonzero determinant. The complex nature of Bij is essential for 
Bij + Bji to be nonsingular. It is important to be cautious when implementing ap- 
proximations for the diffusion matrix and not to prematurely throw out the essential 
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imaginary pieces. Finally, the fact that B12 and B21 are complex implies that these 
cross-correlations cannot be understood on a purely classical basis. 

According to the conventional stochastic interpretation we should view (^) and 
(p7|) as Langevin equations for the classical stochastic variables Xl and tt^ with 
and F| viewed as c/asszca/ noises |jl9|. However, it is not at all obvious why this should 
be true. The Langevin equations ( P^D and ( P7D are operator equations and we must 
have further information about the quantum state of the system before any classical 
interpretation can be accepted. Furthermore, the noise operators do not commute: 



27T 



F^{xr,r],),F^{x2,m) = TT^ks' ks ^^^(r/i - r/2). (37) 



sin ksR 



ksR 



In principle the noises are certainly not classical: we will go on to show that neglect- 
ing the quantum correlations buried in the noises produces results conflicting with 
standard field theory. (The question of the quantum state and the quantum nature 
of the coarse-grained fields has been taken up in more detail in Ref. |^ where it has 
been pointed out that the coarse-graining by itself does not lead to a set of classical 
equations.) Here we follow a different path by deriving a quantum stochastic Liou- 
ville equation that incorporates, at least to some extent, the correlations between the 
noises. 

The fact that F^ and F| do not commute implies that ( [2B| ) and (^) should 
be treated as two separate Langevin equations. Strictly speaking it is not valid 
to substitute (p6D in (^Tl) and treat the resulting equation as a stochastic equation 



second order in time. This will lead to wrong answers for averages involving ttl. The 
approach we will follow avoids this pitfall. 

The Hamiltonian equations of motion (p6|) and (|27D are exact as no approximations 
have been made so far. The first approximation we make is to drop the spatial 
derivative term in (pTf); this is because we will be interested only in the behavior of 
the quantum field at "large" scales. The coarse-graining will be implemented in the 
sense of a temporal ensemble, i.e., we focus attention on one spatially fixed coarse- 
grained domain and consider the evolution of the coarse-grained quantum field defined 
on that domain. The spatial coarse-graining and this interpretation imply that all 
two-point objects be evaluated with the spatial separation between the points being 
much less than the coarse-graining scale, i.e., R -C 27ikg^. With this limit in place 
and with the neglect of spatial derivatives, the quantum Langevin equations are 

(38) 
(39) 




where 



and 



^^(r^i)^- (^2)) ^ 25,,(r/05(r7i - m), (40) 



Bnim) ^ ^fc/ ks \X,,{m)\\ (41) 



i?i2(r7i) ^ ^5 Xfe,(r7i)X,^(r/i), (42) 

B2i{Vi) ^ ^ks'\ks\Xks{Vi)Xl,{vi), (43) 
i?22(r/i) ^ ^fc/|A;5||Xfc,(r7i)f . (44) 

Spatial variations within one coarse-grained domain cannot be sampled by the coarse- 
grained field; this accounts for the fact that there are no terms refiecting such a 
dependence in (|8D-(|||). 

The dynamical equations (|38[) and (^) can just as well be obtained from the 
stochastic Hamiltonian 

H{Xl, ttl) = ^-Kj^ + ia;'(r/)x/ + F^ttl - F^Xl, (45) 

where the time dependent "frequency," 

00^ = S^m^ - |. (46) 

The coarse-grained field is now viewed as the coordinate variable in the one-dimensional 
quantum mechanical problem specified by (^51) . The terms containing and F2 are 
taken to represent stochastic external perturbations with correlations specified by 
(^)-(|4^). (The Hamiltonian (^5]) is a time dependent generalization of the randomly 
forced oscillator considered previously in a different context by Merzbacher pO|.) The 
idea now is to study the one-dimensional quantum mechanical problem instead of the 
original field theory. It is important to note that for a quantum analysis we cannot 



just use the equations of motion (|38D and (|39|) ; a Hamiltonian is necessary. On the 
other hand were we only interested in a classical analysis, the equations of motion 
would suffice. (A discussion of this point is given in Ref. [plj] .) 

Before proceeding further some cautionary remarks are in order. First, while the 
above assumption is an improvement on previous work to the extent that we are 
not assuming the system to be classical, it still does not constitute a well controlled 
approximation scheme. In particular, the Hamiltonian (|45|) has been written down 
simply by fiat. Nevertheless, to see whether the results and insights obtained using 
this approach are persuasive, our attitude will be to take the formalism as it stands 
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and proceed as far as possible without any further assumptions. Second, there is a 
coordinate dependence inherent in the phase space formahsm we will be employing 
shortly: the distribution function is not invariant under canonical transformations 
(this feature is generic to quantum mechanics and is not specific to our problem). We 
will return to these problems in more detail later on. 

The quantum system is completely described by its density matrix, which written 
in the coordinate representation, 

P(^L, X'j} = Y: W.^P.iXLWjiX'^), (47) 
j 

obeys the quantum Liouville equation 

zp{XL, X'J = [HiXL) - H*ix'MXL, X'J. (48) 

The passage to a quantum phase space is now made via the Wigner transform p2| of 
the density matrix: 

fw{XL,PL)= f^e'P^^^p{XL + XL/2,XL-XL/2), (49) 



where the new variables 

Xl = (Xl + X'J/2, (50) 
XL = Xl-X'l- (51) 

The Wigner function fwiXiyPi) is always real and properly normalized over phase 
space (for bounded systems), moreover it is square integrable (a property not shared 
in general by classical distribution functions): 

J dX^dpL fwiXL,PL) = 1, (52) 

( dX^dpL f^{XL,PL) < 7^, (53) 

where in the second expression the equality holds for pure states. Quantum expec- 
tation values for functions of Xl and ttl alone are given correctly as phase space 
averages with respect to the Wigner function, as for example, 

{h{lTL)) = J dXLdpL h{pL)fw{XL,PL) (54) 

but not for mixed operators such as X^^ fc^. This is related to the ordering problem 
in quantum mechanics; the Wigner formalism is associated with Weyl's rule for the 



ordering of operators [23|. A further obstacle to the literal interpretation of a Wigner 
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function as a true distribution function over a classical phase space is the fact that 
in general it is not positive definite. (Fortunately we will not encounter the ordering 
problem nor the lack of positivity in our example.) More on the Wigner function can 



be found in the reviews of Hillery et al and Narcowich p5 |. 

The Wigner transform of the quantum Liouville equation (^81) yields 

d 



dfj 



(55) 



where the Liouville operator has been written as the sum of a systematic piece 



Lo = PL 



d 



dXj 



_d_ 

dpi 



(56) 



and a stochastic piece 



dXi 



dpL 



Xl, Z2 = Pl)- 



(57) 
(58) 



We now implement the strategy of Kubo in order to obtain a simple derivation 
of the stochastic quantum Liouville equation (cases more complicated than the one 
considered here are treated elsewhere |21|). To begin, we focus attention on the 
dynamical effect of Ls by shifting to the interaction picture: 



/^(Xi,pL;r/)=e-^°V(Xi,pL;^). 



(59) 



In terms of a{XL,pi;ri) the Liouville equation ( pSf ) becomes 



d 

—aiXL,PL;r]) 



-e^°''Lse-^°V(Xi,Pi;r/) 
n{r])a{XL,PL;v)- 



This equation has the formal time ordered exponential solution 

a{XL,PL;v) = 



1+ dr]i n{r]i) 

•'VO 

exp ( / drj' ^{rj'] 



V rvi 
drji / dri2 Vt{r]i)Vt{ri2) + 



cr{XL,PL]Vo), 



(60) 
(61) 



a{XL,PL;Vo) 
(62) 



T 



where the initial value a{XL,pL; rjo) is specified at some initial time rjo. All the noise 
terms come multiplied together in each term of the series. If we take the average 
over noise of (p^), these terms will either be zero, or will produce delta functions. It 
is easy to see that only the quadratic product of noise terms needs to be computed, 
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this following from the Gaussian nature of the noises. With ( )^ denoting an average 
over noise, we find, 



dzjdzi 



-Lorn 



(63) 



The noise averaged version of the time ordered exponential solution ( |62D then turns 
out to be 



rv 

{a{XL,PL;v))N = cr{XL,PL;Vo)+ / dr]i B^j{r]i 



dzidzj 



-Lorn 



{a{XL,PL;Vi))N ^ 
(64) 



which may be immediately differentiated to yield 



d 

— {a{XL,PL;v))N = Bijiv) 



dzidzj 



-Lov 



{a{XL,PL;v))N- 



(65) 



We recall that the transformation to the interaction picture involved only Lq which 
is of course unaffected by averages over the noise. Therefore, there is no difficulty in 
writing ( |65D in terms of the original distribution function: 

{fw{XL,PL]V))N = {fw{XL,pL;r]))^ + Bij{r])-^-^ {fw{XL,PL]V))N ■ 

(66) 

This is the required stochastic quantum Liouville equation and, as is obvious, it 
has the standard Fokker-Planck form. Alternatively, ( |66D may be written in a more 
convenient form in terms of the explicitly symmetrized diffusion matrix Dij = {Bij + 



d_ 
di] 



d 1 
— {fw{XL,PL;v))N = -Lo {fw{XL,PL]v))N + 2^'^^^^ dz dz ■ 



{fw{XL,PL;v))N- 
(67) 

The stochastic equation ( [B7| ) and the nature of its derivation merit a few clarification- 
ary remarks. First, it is not necessary to begin with the Wigner formalism; we can 
just as well employ the density matrix (either by following the procedure used here 
or the influence functional approach ||2^). The stochastic equation for the density 
matrix can then be converted to one for the Wigner function by implementing the 
"twisted product" [^. Second, as only the case of free fields is treated here, the 
Hamiltonian is at most quadratic in the dynamical variables. This is why ( |67[ ) is of 
the standard classical Fokker-Planck form; such a simplification does not obtain in 

is classical, this does not imply that 



general [^]. Of course, even if the form of 
all solutions be classical distribution functions. It should be emphasized though that 
Xl and pl are not operators and can be treated as ordinary classical objects. 
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We point out that there is no need to invoke any ad hoc thermodynamic analogy 
(e.g., fluctuation-dissipation relations) in our derivation of the stochastic quantum 



Liouville equation as was done by Graziani in a first attempt to apply the Wigner 
function formalism to stochastic inflation. As we will show in the following sections, 
such relations do not hold in general and any analogy with conventional Brownian 
motion must be treated with extreme caution. A related remark is that since ( p7D is 
formally a master equation one might expect to define a suitable entropy satisfying 



some variant of the H-theorem p9|. For example, it is easy to see that because of 
the diffusion term, the "linear entropy" or "mixing parameter" Trp^ = / dX^dpi 
will always decrease with time (implying that the quantum state is getting more and 
more mixed). This must not be interpreted in the sense of "quantum decoherence" 
pO| as we are dealing with a /ree theory and there is no coupling to some external 
environment. (One way to understand this result may be that this decrease simply 
mirrors the loss of information inherent in our time dependent coarse-graining.) 

Finally we draw attention to some technical issues. Note that no assumption is 
needed as to the symmetry properties of Bijij])] this allows for the fact that the noises 
do not commute. Note also that while separately -B12 and -B21 need not be real, they 
appear in ( p^ ) only in the symmetrized combination 512+-B21 (since partial derivatives 



commute), which, as is clear from (|^) and (43) is always real. The derivation of 
is also free from any kind of "slow-roll" assumption although this merit is mainly 
technical as physical results at late times remain unaffected when such conditions are 
imposed (see Ref. [|], Appendix A). 

We now turn to the problem of solving the stochastic quantum Liouville equation. 
Formally, the solutions are not difficult to obtain as ( |67D is just a Kramers equation 



describing a time dependent Ornstein-Uhlenbeck process The average values 

satisfy 

-^k^i) N = ^iM) ^h) N (68) 

given the initial condition {zi{r]^) ^ = Zio. The matrix Aij is defined by 

d 

Lofw = "^ijg^, i^jfw) (69) 

and in our case. An = A22 = 0, A12 = 1, A21 = —^"^{7]). The propagator for the 
average values (z,), Gij, satisfies, 

^Cjj = AikGkj', Gijirio) = 6ij. (70) 

The second moments follow from 

^ {z,z,)^ = A,, {z,z,)^ + A,, {z,z,)^ + A,. (71) 
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This equation is also obeyed by the covariances Eij = {ziZj) ^ — {zi) ^ i^j) m' which 
can themselves be expressed in terms of the propagator as 

S(r/) = G(r/) ■ S(r/o) ■ G(r/) + /' dr^' G(r/) ■ ^-^(r/') ■ D{r^') ■ ^-^(r/') ■ G(r/). (72) 

The first term is the systematic contribution arising from a reversible dynamical 
evolution from the given initial condition. The second term represents the irreversible 
stochastic contribution due to the diffusion matrix. 

It can be shown ||31| that the general solution of (|67|) with delta function initial 



conditions W{z, rjo) = Hi ^(-^j — -Zio) is 

Wiz,v;z,,vo) = ^{Det S)-i/2exp (-1(5 - (z)^) ■ S'^ ■ (z - (^)^)) . (73) 

The function W{z, r/; zq, rjo) serves as the propagator for the Fokker-Planck equation 
(|67|). Solutions for arbitrary initial conditions fwi^o^Vo) can be generated from it by 

fwiz, V) = J dzo W{z, T]; zq, Vo)fwizo, Vo)- (74) 

These general results will be applied to exponential and power law expansions in the 
following sections. 
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III. Exponential Inflation 

The case of a de Sitter expansion furnishes a particularly simple example in which 
to implement the procedures of Sec. 11. Our stochastic approach not only reproduces 
some previous field theoretic results but also introduces a new interpretive framework. 
In this section we aim mainly to obtain quantitative results. 

The scale factor for de Sitter space is 

SM ^ (75) 

and the conformal time 

V = (76) 

We note that here we are not really interested in the case of an eternal de Sitter 
expansion. Initial conditions for stochastic inflation will be assigned in the finite 
past, at the beginning of the inflationary phase. 
The scalar field modes now satisfy 

.2 \ \ 

0. (77) 




If the mass is zero or at least small compared to Hq^, the "unstable" sector is char- 
acterized by A;^ < 2/?]^. Therefore, following Starobinsky |]T|] it makes sense to set 

ksiv) = (78) 

where e is a constant that serves to parametrize the cutofl. If we assume that inflation 
began at the time r]o, this implying a natural infrared cutoff r]Q^ (more details may be 
found in the Appendix), it is clear that e cannot be arbitrarily small as we must have 
g^-i ^ ^-1 jf Qj^g jg interested only in late time results, i.e., when rj <^riQ then e may 
be taken to be small. In this paper we will not restrict e to be arbitrarily small but 
will allow it to be as large as unity. In principle, it is desirable that physical answers 
not depend on e; this will turn out not to be the case. As will be shown later all 
infrared divergent quantities are only weakly dependent on e but this does not hold 
in general (the situation is more complicated for power law inflation). Our stochastic 
approach will correctly reproduce the cutoff dependences for infrared flnite quantities 
calculated from conventional quantum fleld theory (see the Appendix). 

A curious special property of de Sitter space is that even when the mass is non- 
zero (no infrared divergence), there is still an initial growth of ($^) to an asymptotic 
limit ($^)^^, which is the value in the Bunch-Davies vacuum [|15]. The mathematical 
reason for this behavior is simply that the mass and curvature contributions in ( [77D 
scale identically with conformal time and that for m? small compared to there is 
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still an "unstable" infrared regime despite there being no infrared divergence. As will 
be made clear in the next section this feature is not shared by power law inflation. 
The mode equation (^) admits the general solution 



(79) 



where 



9 



4 H^- 



(80) 



The arbitrariness of the de Sitter vacuum is reflected in the various possible choices 
for Ci and C2. In this paper the quantum state we will use is the Bunch-Davies 
vacuum |T^, characterized by Ci = and C2 = 



The symmetrized diffusion matrix Dij now follows from (p|)-(]4^): 



OTC L 

- D2i{r]i), 

/I 

D22(r/i) = 2522(r7i) - —Vi' 



(81) 



in 



^2) 
^3) 



Our eventual goal is to solve the Fokker-Planck equation (^) given the diffusion 



coefficients (Bl)-(83). The potential term in the systematic component ( p6D of the 
stochastic Liouville operator is characterized by 



u 

4 



1]^ 



^4) 



For all the cases we consider in this paper uj'^{t]) will be negative (as z/^ > 1/4). We 
are dealing therefore with a time dependent upside down harmonic oscillator. The 
equation for the propagator ([T0|) can now be easily solved, and in terms of the initial 
time Tjo, we find 



G22{V) 



1 



1^+1/2 



-u+1/2 



+ {u + l/2) 



a/2 



2u 



{u^ - 1/4) 



2z/ 



1 

2^ 



{mo)-"' 



V 



.10) 

Vo, 

u-l/2 



(85) 
(86) 
(87) 
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This is the general solution, valid for all the special cases considered in this paper. 
We now confine attention to the massless case where the parameter v = 3/2, and 



71 I krj 



r(2)/; N _ l'^^^^ 



1/2 —ikr) 



H^Mkv) = (^J J^-.i^-^v). (90) 



The diffusion matrix then takes the simple form, 

(^i) ^ ^^r'(i + e'), (91) 

D22{vi) ^ i^^r' (l - + • (93) 

Note that to leading order (for e <^ 1) the diffusion matrix is independent of e. It 
is misleading however to conclude that the actual value of e is not important as at 
this order Det Sij = (notice that this is a direct consequence of retaining noise 
cross-correlations). In order to eventually obtain a nonsingular covariance matrix we 
must go beyond this level of approximation; the final solution for the Wigner function 
is in fact strongly dependent on e. 

It is a tedious but straightforward exercise to obtain the covariance matrix using 
(^). We assume an initial distribution such that Zio = but impose no conditions 
on Sijlrjo). Ignoring for the moment the systematic component, the stochastic con- 
tribution turns out to be 



187r2r;o 



^22 





g2 ^2^3' 



3 3r7o' 



(95) 



g2 g2^3^ 



(96) 



One can easily check using (^) and (|5BD-(P5D that at late times the systematic 
contribution to the covariance matrix is negligible compared to the stochastic piece. 
(Initial conditions are discussed further below.) 

The full solution for the noise averaged Wigner distribution function follows triv- 
ially from ( [73D and (0). In this section we will concentrate only on the covariance 
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matrix itself as all average values of interest can be computed directly from it. De- 
tailed study of the distribution function will be postponed to Sec. V. 

The covariance matrix (p^)-(|96D refers to the "conformal" variables and p^. 
Reverting to the original field $, we introduce new c- number variables 0c and Pc via 
the canonical transformation 

0c = (97) 
p, = Spl-SXl. (98) 

The corresponding covariance matrix may be written as 

(0) _ /j. 2\ _ c-2 



^11 



= S-^Eu, (99) 



^12 - -21 - {<PcPc) N = -12-^2ii, (100) 



E'^2' = {p/)j, = S'E,2~2SSE,2 + S'Eu, (101) 

where all averages of the type {ziirj)) vanish as a consequence of our choice ZiQ = 
for the initial condition. Comparison with the field theoretic results is simpler if 
we introduce the "velocity" 0^ (the prime denotes differentiation with respect to 
cosmic time) in place of the canonical momentum pc- Noting that for de Sitter space, 
In-qo/t] = HQ{t — to), where to denotes the beginning of the exponential expansion, 
and using (p4D-(p^), the new covariance matrix (for a massless field) turns out to be 

^2\ 



^cYc/N 




e 



(102) 



3 6 

= If^S^-e-'"""-""). (104) 



£1 ^ Le-3H.(<-<.) I _ (103) 



The above results record only the stochastic contribution to the covariance matrix. 
It is easy to compute the systematic contribution for an arbitrary initial choice of 
Eij from ( [72D and (|85D-(p8D (since the Wigner distribution function must be square 
integrable we cannot take the initial distribution to be a delta function over phase 
space). The contribution to Eu consists of a constant piece and terms that fall off 
exponentially with cosmic time. Contributions to S12 and E22 also display a similar 
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exponential fall-off. It is important to note that these contributions, though insignif- 
icant at late times, can dominate similar terms that already exist in the stochastic 
piece (especially for small values of e). Therefore, only the late time limit is inde- 
pendent of initial conditions. This is in contrast to the field theoretic case where the 
contribution from initial conditions is usually irrelevant even at early times (see the 
Appendix). We also draw attention to the fact that the exponential fall-offs in the 
stochastic calculation are not the same as the field theoretic ones; again, this is of no 
consequence at late times. 

As long as the initial distribution is such that (0c) at ^"^^ {Pc)n are zero (i.e., 
^io = 0), at late times [i] small) and with e <^ 1, the leading order contributions are 



"11 - \Vc 



=2^ = /0'2\ ^ ^ -"O _ HOT) 
S6 24772 y^"") 

The E^i term reproduces the standard quantum field theoretic result |JTB|] (and ( p.66| ), 
the Appendix) for the expectation value ($^), here viewed as a noise average for the 
c-number variable 0c provided that e is small. However, if we set e ~ 1 the answer 
does not agree with the field theoretic result ( |166|) found in the Appendix (which 
unlike the stochastic calculation is essentially cutoff independent). 

For the sake of comparison, if we set e = 1 in (|102|) - (|104|) we find, at late times, 



N 367r2 



■^11 



it -to), (108) 



KIT 3 

^ = (<Pc€).-^^, (109) 

TIjW TT 4 

^22 _ / ,/ 2\ ^ -"O 



S' V- % - 24^5- 

whereas the corresponding field theoretic results ( |166| ), ( p.70| ), and ( p.73| ) of the Ap- 
pendix give, at late times: 

^ %it-to), (111) 



47r2 

1 f2rr 3 

_/$$' + $'$) ~ lii^ (112) 

2 ^ ^ 87r2 ' ^ ' 



~ lil^. (113) 



^4 rr 4 

167r2 

We see that while the stochastic results for and correctly reproduce the 
the cutoff dependence found in the field theoretic case, the numerical values of the 
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coefficients do not match. This is not a serious problem as these quantities need to 
be renormahzed anyway (something that is beyond the scope of this paper). 

The technical reason for the disagreement between the field theoretic and stochas- 
tic calculations is the neglect of spatial derivatives in (pOf). One is attempting to 
approximate a time dependent quantum sum over modes by a modified dynamics 
(via the noise term) for the zero mode and neglecting all the other modes (apart 
from their contribution to the noise). When computing ($^) via the standard field 
theoretic method the infrared sector provides the dominant contribution. On the 
other hand, when computing ($$' + $'$) /2 and (<l>'^) extra multiplicative factors of 
k and fc^ (see the Appendix) weaken this infrared dependence. One expects therefore 
that the stochastic method should work better for (<l>^) and, as we have seen, this is 
indeed the case. The fact that stochastic results are more accurate for small values 
of e is also easy to appreciate: a small e means that only long wavelength modes are 
contributing to the noise so that k is indeed small and can be neglected. However, if 
we let e be of order unity, then the neglect of spatial derivatives will lead to errors. 
A possible remedy is to work with a Wigner functional defined directly from the field 
theory but this may well be at the expense of the calculational ease that characterizes 
the present approach. 

We now consider the massive field but confine attention to the case <^ Hq^. 
The propagator G{ri) is still given by (|85|)-(|88D except that the parameter u is now 
given by 

(114) 



3 
2 



m 



In the limit of a small mass the diffusion matrix is essentially the same as for the 
massless case (since z/ ~ 3/2). Keeping the diffusion matrix given by (pTD-(p3D 
but using (|114|) for the propagator it is a simple matter to solve for the stochastic 
contribution to the covariance matrix. The final expressions are very long and not 
very illuminating. Here we present only the leading order stochastic terms at late 



times: 



" /N 



N 



0. 




(115) 
(116) 



;il7) 



where we have dropped all terms that vanish faster at late times and also neglected 
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terms that are of higher order in /Hq^ . An interesting feature of the massive case is 
that all contributions arising from ^•^■'(r^o) are negligible even at early times (unlike 
the massless case). The agreement with field theory is remarkably good. With e 
small, or more precisely, in the range. 



exp(-3ifoV(2m2)) < < / H^^ 



(118) 



all asymptotic late time values are exactly reproduced (similar inequalities are derived 
somewhat differently in Refs. and |[T^) . Unlike the massless case, this time the 
approach to these late time values is also in agreement with the field theoretic results 
(i.e., no mismatch in the exponentially falling off terms). At late times the above 
expressions reduce to 



"11 

"12 

"22 
^6 



" /N 



kl 2 



N 



N 




(119) 
(120) 
(121) 



in agreement with the field theoretic results (|178|) , (|182|) , and ( |184| ) of the Appendix 
(with e in the previously indicated range). Unlike the massless case there is no leading 
order e dependence in S12 and H22. To avoid Det S,-^^ = 0, it is important to keep the 
sub dominant mass squared term in ( |119| ). It is again easy to verify that agreement 
with field theoretic results does not extend to the case e ~ 1; the two calculations 
now differ by multiplicative factors of order unity. 
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IV. Power Law Inflation 

In this section we treat a power law expansion a{t) ~ with p > 1 and consider 
only the case of a massless scalar field. Assuming that infiation set in at the initial 
time rjQ (with the scale factor set to unity at this time), we find from (0), 

/ x(l-2.)/2 

Siv)= -] (122) 



where 



It then follows that 



1 -3p 

2(1 -P)' 



(123) 



It is useful to note that for p > 1, z/ > 3/2 and also that as t ^ oo, rj ^ 0. In the 
formal limit p oo, u = 3/2 which is the value for de Sitter space. 
The scalar field modes now satisfy 

Xfc + I F + ^ 1 +[l-u']^\Xk = 0. (125) 




This time the mass and curvature contributions scale differently with conformal time. 
It is easy to see that at late times the mass term dominates the curvature contribution 
since u > 3/2. Therefore, as cosmic time increases (and the conformal time decreases) 
there is a continuous fiow from the "unstable" to the "stable" sector. However, 
destabilization will still occur for the massless case where the "unstable" sector is 
characterized by k"^ < (z/^ — l/A)/iff. We see also that even in the case of power law 
infiation it is only natural to implement the same choice that we made in the last 
section, i.e., to set ksijf) = e/rj. 

The mode equation ( |125|) for a massless field admits the general solution 

Xkiv) = C,r,''^H'^\kr,) + C^ri^'^H^Kkr,) (126) 

or, in terms of the original field, 

UV) = C^r,,"' H^^\kv) + C.r^o'" Hl'\kv). (127) 

The adiabatic vacuum is specified by Ci = 0, C2 = \^/2, i.e., 

1 /2 

MV) = (y) Hi'\kv). (128) 

Assuming that all "high frequency" modes were in the adiabatic vacuum at the onset 
of infiation, ( |128| ) enables us to compute the diffusion coefficients (|^)-(PD which are 
the same as (|8T|)-(|83D except that now u is specified by ( |123|) . 
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Power law inflation with p ^ 1 can be treated in a simple and direct manner 
by following the same approach as that for the massive field in exponential inflation. 
Note that when p is large, 



3 111 



(129) 

p p^ p-' 

This allows us to approximate the diffusion matrix by the one for v = 3/2, (^T])-(53). 
The propagator G{ri) is given by (|85|)-(|88D with u specified by ( |129|) . The stochastic 
piece of the covariance matrix can now be found by a straightforward computation. 
The result is too long to write out in entirety and we content ourselves by just 
displaying the leading order terms: 



--11 



^12 



" In 



S3 



cl N 



^22 



2 




(130) 
(131) 
(132) 
(133) 

(134) 
(135) 



where terms lower order in 1/p and vanishing faster at late times have been dropped. 
The variables (pc and Pc are still defined by (^) and ( P5| ) except that S{ri) is now given 
by ( |122|) . For small e, (|131|) is in agreement with the field theoretic calculation of [|16 



(also compare with ( |190| ) of the Appendix). However, just as for the massive field 
in de Sitter space, e cannot be arbitrarily small. Consistency with the field theoretic 
results (|1971) and (|203|) of the Appendix requires that 



p ^ < < 1 



(136) 



in ( |133| ) and (|135|) . The above expectation values were also computed by Kandrup 
using a different method 0. While our result for (0^^) is in agreement with his, 
this is not true for the other two cases. The inconsistency can be traced to an 
approximation for the noise that does not take the commutator properly into account 
(see the discussion of this point in Sec. II). 

In this case, while at late times (0c^)iv go^s to a constant (as for the massive field 
in de Sitter space), {(Pc(Pc)n {4>'c'^)n vanish. This is in contrast with the case of 
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a massive field in de Sitter space where these quantities instead of vanishing, also go 
to constant values. The role of initial conditions is similar to that for the massive 
field in de Sitter space rather than the massless one: the systematic contribution 
to the covariance matrix is always negligible as long as p ^ 1 (but not otherwise). 
Consequently, the late time results follow from ( |131| )- (|135|) : 



- £^^o-^ (137) 



^2 



to 



{<f>c€)N ^ Y2^2^o'[j) , (138) 

where we have taken e <^ 1. For e = 1, we have, 

;Vo"^ (140) 



2\ 13P _2 



' 'N 727r2 

(M'c)n - 4^2^0'{^j) , (141) 



ki 2\ ^. 1 „-4 [to 



> - 24.^"" If j ■ <"^> 

The field theoretic results (|19CI|) , (|198|) , and ( p04|) of the Appendix yield the cor- 
responding late time limits. 



2\ P ^-2 

i(44.' + t't> ^ ^%-'(^) , (144) 



;*^) - i^h', (143) 

1 /^.T./ , ;f,t;f,\ - . ^ -3 fto^ 



to 



There is reasonable agreement with the stochastic results when e is small but as is 
expected the results diverge from each other when e ~ 1. 



The limit p — > oo may be applied to (|131| ) using 



'^V'' = l + ^lnW+..., (146) 
,^0/ P \VoJ 



with the result 

1 9 1 
1 + -e^ + -e 
3 9 



(147) 



2\ ^ Hpit — tp) 
= /n - 4TT^r]Q^ 

Since Hq = 1]^^, this agrees with the result ( |102| ) for a massless field in de Sitter 
space. In a similar manner one can check that {4>c<P'c)n Wc"^) n reduce to the 
appropriate expressions for a massless field in de Sitter space as calculated from the 
stochastic approach. 
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V. Solutions and Interpretations 

In this section we study tlie full phase space distribution function. Given that 
we have already computed the relevant covariance matrices it is now a simple mat- 
ter to write out the corresponding Wigner functions. In the examples studied here 
these distributions will be positive definite and as such may be interpreted as true 
probability distributions, at least formally. 

A knowledge of the distribution function is important as it will enable us to 
critically address issues such as the existence of fluctuation-dissipation relations and 
whether there exist late time thermal solutions or not. These are the problems we 
will tackle first. 

The stochastic Liouville equation ( |67|) is written in terms of the conformal variables 
Xl and p^. In all the examples we studied, = 0, in which case the general solution 
([f3|) becomes, at late times, 

fw{z, v) = ^{Det exp (^-^~z ■ E'' ■ , (148) 

all contributions from nontrivial initial conditions having washed out in this limit. 
Converting to the variables (pc and Pc appropriate to the original frame, the distribu- 
tion function ( [I48| ) goes over to 



mc,Pc) = ^iDet sW)-V2exp (-]- 

ZTT \ Z 



"11 9c ^ ^"12 VcPc ^ ^22 Pc 



(149) 

It is important to appreciate that while fd gives the correct expectation values (Sees. 
Ill and IV) and is a perfectly respectable classical distribution, it is not a Wigner 
function defined from the beginning for the variables 0c and vr^. This is because, as we 
noted earlier, these distributions are not invariant under canonical transformations. 
(We are treating the conformal variables Xl and ttl as the preferred variables to 
quantize.) However, the key point is that in our case the linear entropy remains 
invariant under this transformation. 

A key observation regarding ( p.49| ) is that a knowledge of the reduced distribution 



/r(0c) = / dpcfw{4>c,Pc 



00 



1 



^11 




(150) 



is of no use in reconstructing the original distribution. This trivial fact has important 
consequences if one attempts thermodynamic interpretations of the results from our 
stochastic analysis using only the reduced distribution. Other points to keep in mind 
are that, in some cases, to leading order in e, Det S^"^) = (therefore the distribution 
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function is not independent of the cutoff), and that the cross-term proportional to (pcPc 
represents a nontrivial contribution from quantum correlations. In order to discuss 
these issues more concretely we now return to the specific cases studied earlier. 

We consider first the massive free scalar field in an exponentially expanding Uni- 
verse. At late times, with e satisfying the condition ( |118| ), we have. 



6 , 1 (9Hn 



' "(151) 

Clearly this distribution is not stationary because of the dependence on Sir]). On the 
other hand, the corresponding reduced distribution has the "equilibrium" form 

E^' V 3 """" V ^ 



Mc) = -T7-2\h^^P\-^nrT^"] (152) 



(2;^) ^''Pi-f^OHEHi'Pc)) (153) 



where Pq}j = Tgh = Hq/2'k is the Gibbons- Hawking temperature of de Sitter space 
||3^ , Vh = AtiHq^ /?) is the three- volume within the Hubble radius, and Eh{(Pc) = 
V{(f)c)VH is the energy of the scalar field within that volume (the kinetic energy is 
not important if the field is in the "slow-roll" regime). The thermodynamic interpre- 
tation of the stochastic formalism P] was suggested by the striking Boltzmann-like 
nature of ( |153| ). However, the full distribution ( |151| ) does not seem to encourage such 
speculation: it is not stationary nor of the form e~^^ (nor is the reduced distribution 
friPc) of the form e~''^'='"). 

We observe that transforming to a new variable Vc = Pc/S^ makes the late time 
distribution ( |153| ) time independent. However this transformation is not canonical 
and does not preserve the linear entropy. Therefore the distribution f{4>c,Vc) is not 
physical. In any case such a trick fails for the case of power law inflation: there the 
phase space distribution cannot be made time independent. 

Turning now to the massless case, at late times. 



fcl{4>c,Pc) 



247r2 / 271 

exp 



2 




cPc ^ QPc 



(154) 

The singular nature of this solution as e — > is apparent (as is the fact that it is 
explicitly time dependent). Note, however, that in this case. 



which is independent of e. The late time reduced distribution ( |155|) is a solution of 
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the diffusion equation 

with D = Hq^/Atc'^. This is suggestive of a (usual, time independent) random walk 
interpretation. However, the time dependence of the terms oc p^^ and oc (pcPc in the 
full distribution are hard to reconcile with this view. 

The case of a massless field in a power law spacetime is treated next. Here the 
late time full and reduced distributions are, respectively: 




As with the massless field result (|154|) , the full distribution is again singular in the 
limit e — > 0. Also the reduced distribution ( |158| ) is independent of e, as in the other 
cases. Unlike the other two cases however, it does not seem to have any "natural" 
interpretation. 

The late time linear entropies cr = / d(j)cdpcfj for the three cases studied above 
are respectively: 

<^ds = ; = — 5 , =, (160) 

^ ^%M^(i-to) \ H,\t - to) 



In all cases a is approximately proportional to S~^. The possible significance of this 
result will be discussed later below. 

It is by now clear that the late time phase space distributions obtained here are 
very difficult to fit into a conventional Brownian motion picture. In fact, this is a 
very obvious point and manifest in our stochastic Hamiltonian (^5]). In standard 
Brownian motion the environment with which the system interacts produces both 
dissipative and diffusive effects. The dissipative effects arise from the back reaction 
of the environment. Such an effect is absent in stochastic Hamiltonians of the type 
(^). In principle, then, there simply cannot be a fluctuation-dissipation theorem of 
the usual sort: this conclusion is manifest in the fact that in no case are our late time 
solutions for the distribution function stationary. (However, this does not mean that 
there cannot be asymptotically constant values for some average quantities.) 
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We recall that the origin of the stochastic noise is simply because the "system 
size" is changing with time and not because of some external interaction. It appears 
that a mistreatment of this key point has led some authors to claim that quantum 
decoherence occurs in this model. That in fact it does not can be explained by the 
following direct argument for which the author is indebted to Juan Pablo Paz. 

The transition from quantum to classical was studied in the context of stochastic 



inflation by Morikawa |Q and Nambu |Q who analyzed the properties of the evolu- 
tion operator for the reduced density matrix of the long wavelength modes. When this 
propagator is written in path integral form the effect of the short wavelength modes 
appears to be contained in a term that is rather similar to the Feynman- Vernon in- 
fluence functional |2^, F = expzF. In ordinary open systems, the presence of an 
imaginary part in the influence action F produces a tendency towards diagonaliza- 
tion of the reduced density matrix in a flxed basis. This is known as decoherence. 
For stochastic inflation, the imaginary part of F was calculated and related to deco- 
herence. However, it is possible to show that this interpretation is not correct and 
that there is no decoherence produced by the coarse-graining of stochastic inflation. 
The basic reason is that the time dependent nature of the coarse-graining prevents us 
from interpreting the influence functional in the usual way. In fact the reduced density 
matrix at a given time can be written as the product I\k<ks{t) Pk- As the number of 
modes present in the system varies with time, the evolution operator J{t, to) has some 
peculiar properties. It can be written as a product of an evolution operator for each 
mode Jfc(t,to) where for k < A;s(to) (modes that were already present in the system 
at if: = to), the Jkit^to) are ordinary unitary operators while for ks(to) < k < ks(t) 
(modes that enter the system between to and t) the evolution operator is simply 
Jk{t,to) = Pk{t)- If one writes these operators in path integral form one realizes that 
there are real exponential terms simply due to the fact that, if the state of the fleld is 
the vacuum, the reduced density matrix pr(0fe, (p'k) is a Gaussian. The only effect that 
the "influence functional" has in this case is to generate the above Gaussian factors. 
It is clear that this is not related to decoherence but to the fact that new modes are 
entering into the system and that the evolution operator fully contains the reduced 
density matrix of the incoming modes. 

Another argument put forward for a late time classical limit was that since the 
commutator (^) is oc e^, it is "small" and can be ignored. This of course cannot 
be correct. The reason is that it is not just a single mode commutator one has to 
look at but the total integrated contribution from the initial time to the flnal time of 
interest. This is not a negligible fraction. 

Is there a classical limit or not intrinsic to the formalism? The linear entropy 
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does decrease exceedingly rapidly as shown by ( |159|) - (|161|) but it is not clear what 
this means: we have just argued against interpreting this sort of decrease as being 
due to quantum decoherence. An intuitive basis for this result may be that it reflects 
the loss of information inherent in our time dependent coarse-graining. With the 
passage of cosmic time, two-point functions are averaged over ever smaller comoving 
volumes. The "smearing" scale is set by kg-^ = v/^ ^ind r] ^ as t ^ oo. Since in 
our formalism we are tracking only one coarse-grained domain throughout its history 
this represents a loss of information with cosmic time. We may speculate plausibly 
that the decrease of a as supports this viewpoint. However, just because our 
knowledge is incomplete is no reason to suppose that the Universe is becoming more 
classical! 

At the present stage of analysis and understanding it appears unlikely that the 
quantum to classical transition in the early Universe can be explained by the stochas- 
tic paradigm. In particular, the treatment of density perturbations by modeling 
quantum fluctuations as classical noise appears to be unjustified. 
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V. Conclusion 

This paper's main concern was to model a free field theory in an inflationary Uni- 
verse by way of a stochastic quantum Hamiltonian. It was shown that this approach 
produced results that agreed well with those from straightforward quantum field the- 
ory. Furthermore, the role of the length scale parameter e and of initial conditions 
was considered more fully than in previous work. 

The quantum phase space distribution used in this paper enables a consideration 
of quantum correlations that would otherwise be missed. As a result we find that the 
e — limit is singular as far as the distribution function is concerned. This means 
that a finite value of e is essential for the formalism to make sense and that contrary 
to previous belief this parameter does not drop out of the problem. The full phase 
space distribution also enables a critical assessment of such issues as the existence of 
fluctuation-dissipation relations. We showed that fluctuation-dissipation relations do 
not hold (as indicated by the fact that the late time solutions are not thermal or even 
stationary). However, at least in de Sitter space, the reduced distributions for the fleld 
variable alone have very suggestive forms corresponding as they do to a Boltzmann 
distribution at the Gibbons-Hawking temperature for the massive fleld, and to a 
"random walk" distribution for the massless fleld. No such simple distribution appears 
in the case of a power law inflation. The signiflcance of these results remains unclear 
at present. 

We found that in order to obtain results more or less consistent with conventional 
field theoretic calculations quantum correlations could not be neglected. It was also 
pointed out that quantum decoherence does not occur in the stochastic approach. As 
a consequence of these two results, the quantum to classical transition in the early 
Universe does not seem to be intrinsic to the stochastic approach. Directly modeling 
quantum fluctuations by classical noises as a way to study density perturbations from 
infiation is therefore a questionable enterprise. 

There are of course many unanswered questions, chief among them is what hap- 
pens when interacting fields are considered and back reaction is included. This we 
leave to future work. Furthermore, while it is true that the stochastic model "works," 
at least to some extent, we have stressed that it is not free from interpretational 
problems. One can only speculate whether insights gained from this approach will 
actually turn out to be valuable when the full quantum field theoretic computations 
are eventually done. 
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Appendix 

A brief review of conventional field theoretic computations of the various expec- 
tation values of interest will now be given. This will enable us to check results from 
the stochastic analysis. To obtain finite results we will impose an infrared cutoff in 
momentum space at k = tjq^ and an ultraviolet cutoff at k = erj^^. A discussion of 
the reasons for picking these cutoffs will be given at the end of the Appendix. 

We begin with exponential inflation and treat the two cases of a massless field 
and of a massive field with a small mass {m? <^ -f^o^)- Considering first the case of a 
massless field, the parameter z/ = 3/2 (from (^)), and for the Bunch-Davies vacuum 
the mode functions (j)k = Xk/S are 

Uri)={^j)'' H,r,^''Hf/,{k^). (162) 
Consider first, the equal time expectation value. 



{mm) = ^ J^_, dk k'^^ \<PM\ , (163) 



where R= \x — y\. Using the exact form of the Hankel function (|90D it is easy to see 
that 2 

10.(^)1' = §j[l + W], (164) 

hence the integral in (|163|) is infrared divergent and an infrared cutoff is necessary. 
With < 1, we find 

2" 



$2\ = :^ 

/ 4vr2 



(165) 



Notice that since erj ^ > % \ ($^) as computed above is strictly positive. To write 
the result in terms of the cosmic time, we note that tjqI]^^ = exp HQ{t — to), in which 
case 3 2 9 

(*^> = - «o) + ^ ('- + r1 - ^e--.<-«, (166) 

The last term vanishes at late times and the second term is an irrelevant constant 



absorbed in the infrared cutoff. The first term gives the usual result ||I5[ . Notice that 
this term is independent of e; any potential dependence is lost in the infrared cutoff. 
Now we turn to the quantity 

1 / . . \ 1 r''" 



+ $$) = -L 1'"^ dk k'Re {Mvmv)) (167) 



2 

where we have already set kR <^ 1. From ( |162| ) 



4 



Hokr]^/'Hi%{kr]), (168) 
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and it is easy to compute that 

1 



El 

87r2 



(169) 



If we return to the cosmic time, then with a prime denoting differentiation with 
respect to t, 



1 



TT 3 



(170) 



2 ■ ■ svr- 

The last term vanishes at late times and the first term is a constant that is strongly 
dependent on the upper cutoff. (Even if e <^ 1, consistency requires ^ ijiNof 
the first term always dominates.) The late time answer being a strong function of 
the cutoff is simply a consequence of the extra multiplicative factor of k in 0^ which 
not only renders the integral in ( |167| ) infrared finite but also shifts the dominant 
contribution from the integrand towards the upper cutoff. 
To compute i^"^) we first use ( [L68| ) to show that 



(171) 



With kR <^ 1, it is now easy to find 



or, in terms of the cosmic time, 

■$'2 



27r2 

Hi 
167r2 



HI 
167r2 



dk k^ 



e - — 



1 

^0. 



(172) 



-iHoit-to) 



(173) 



The late time value is again a strongly cutoff dependent constant. 

Similar calculations will now be performed for the case of a small mass, i.e., for 
<^ Hq^. Disregarding an irrelevant phase term for real z/, the mode functions are 



now 



TC 



1/2 



Hov'/'Hl'\kr^) 



where, for a small mass. 



m 



z/ ~ 



(174) 



(175) 



2 3i/o^ ' 

Since kr] < 1 over the range of integration t/q^ < k < er]^^, we may approximate the 
Hankel functions by 



Hl'\z) 
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(176) 



and using ( |163|) compute, 

$2^ 



exp 



2m2 
3^ 



As long as e ^ exp(— 3-?/g^/2m^) (i.e., e cannot be arbitrarily small), 

2m2 



1 + 



Ine + 



and at late times 



(177) 
(178) 

(179) 
(180) 



a well known result 15 



The calculational strategy used above can also be implemented to find that 



1 H ''^ 



2m2/3H(,2 



1+ - 



;i8i) 



and that. 



1 



Ha 



^2m2/3i?o^ 



+ e — exp 



2m^ 



{t - to) j (l + e- 



2H„(t-to) 



( 



^2) 



Setting = in (|182|) we recover the previous result 
field. The late time limit for e ^ exp(— 3iJQ^/2m^) is 



2 ^ ^ 87r2 



1 + e^ 



for a massless scalar 



;i83) 



Unlike the massless case, here the result is essentially e independent provided e is 
small compared to one. However, unlike the situation for ($^), the late time value 
does depend on whether e is small or of order unity. 

The expectation value ($'^) can be found in exactly the same way: 



/2 



Hi 

47r2 

4{ 



m 



e^-e 



exp 



2m2 
3^ 



(t - to) 



-4Ho{t-to)\ I 2"^^ r^2 



-2//o(t-to)~ 



;i84) 



If we assume e to be much smaller than unity, then at late times we obtain to leading 
order an e independent result. 



/2 



247r2 ■ 



(185) 
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On the other hand if we set e to be one, then the late time hmit becomes 



12 



167r2 



;i86) 



in 



which is completely different from ( |185|) . It is easy to verify that setting 
( |184|) reproduces the answer (|173|) for the massless case. Therefore for e ~ 1, the 
massless and massive cases give the same result (but not when e is small). 

We turn now to power law inflation and consider the case of a massless field. 
Recall that the parameter u is now given by 

l-3p 
2(1 -p) 

and that the adiabatic vacuum modes are 



(187) 



The expectation value {^'^) is still given by (|163| ) and we can still use the approxima- 
tion ( |176| ) for the Hankel function. If we set = 3/2 then the de Sitter results for a 
massless field are recovered. For z/ 7^ 3/2, we have 

22i/-3 



2\-2v 



7r3(2z/ - 3) 

22--3rrz/)2 



dk fc^d-'^). 

. V J 



.2/{l~P) 



(189) 



7r3(2z/-3) 

reproducing the result of [0. (As long as the power law p > 1, the parameter 
u > 3/2, and it follows that the integral in ( |189| ) is infrared divergent. The lower 
cutoff is necessary to get a finite answer.) We see that ($^) starting from some initial 
value rises to a constant. When the power law p ^ 1, ( p. 8 7]) implies that z/ ~ 3/2 + 1/p, 
and 

(190) 



P -2 



-2/p 



If e >> e"P/^, then 



-2/p 



P 



Ine + 



(191) 



and (|190D is essentially cutoff independent not only as to the late time constant value 
but also as to how this value is approached. The situation is different when p is not 
large. For example, if we consider p = 2 corresponding to z/ = 5/2, 



9 
8^ 



Vo 



-2 



(192) 
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This time while the asymptotic value of (<l'^) is indeed cutoff independent, the ap- 
proach to it is a strong function of e. 

We proceed now to evaluate + /2. The usual procedure, beginning from 
O, yields 



2 ^ 



T{v)T{v - 1)7/0^^2^7/ / _ dkk 



TV- 



4-2u 



(193) 



The integral in ( |193| ) is infrared divergent for u > 5/2 (i.e., 1 < p < 2). For u = 5/2, 



2 



47r2 



-4 1 / ^VO 

Vo ^In I — 



(194) 



which vanishes at late times. In terms of the cosmic time. 



2 Air- 



t 



(195) 



After an initial period of growth, at late times this expectation value vanishes. 
When p 2, we obtain from ( |193D 

t 



1 2^'^""'^ 

2 ^ ^ 87r3(5 - 2u) 



3 / '^0 



(196) 



It is trivial to check that this quantity is always positive. At late times it vanishes 
but the dependence on e is a function of the power law. For large p, (|196|) becomes 



-3 / \ 3/p 

87r2 Uo / 



2-2/p- 



-3 / ^0 



(197) 



:i98) 



If p oo, ( p.97| ) reproduces the de Sitter result (|170|) for a massless field (with 

^0 = ^0^'). 

Finally we compute \^'^)- The standard calculation yields 



^22('^-i)r(z/ - 1) V 



(iA; A;' 



6~2u 



(199) 



The integral is infrared divergent provided z/ > 7/2 (i.e., 1 < p < 3/2). For the 
special case u = 7/2 corresponding to p = 3/2, we have 



or, 



9 -6^2 / £7/0 
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(200) 



(201) 
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At late times this expectation value vanishes. 
In the general case [u ^ 7/2) we find 



22('-i)77o-" 
87r3(7- 2z/) 



-T(u-1] 



+ \ 4 
If 



7-2u 



tQ^2{2p-3)' 



(202) 



which also vanishes at late times. When p ^ 1, ( |202| ) gives to leading order, 
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167r2 yrio 
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r2 



4/p 
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+ \ 4 



(203) 



(204) 



The limit p ^ oo taken in (|203|) gives back the de Sitter result ( |173| ) for a massless 
field. 

We now explain the origin of the cutoffs in the momentum integrals. To prevent 



infrared divergences we follow the strategy of Ford and Parker jl^ by assuming that 
for r]o > V the Universe is radiation dominated and that the quantum state is the 
conformal vacuum. Matching the field modes and their time derivatives at rj = rjo, 
one finds 



1+ £ ikmr'-Tiuy 



(205) 



Since the upper cutoff forces krj < 1, we can use the small argument form ( |176[ ) of 
the Hankel function and compute 



|2 ^ 



^Tiu)Yo-"'k-'nC,ik)~C,ik)\\ 



(206) 



We are now in a position to compute the pre-inflationary contributions to the various 
expectation values of interest. First, consider 
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dk k' 



dk k 







which is a constant independent of u. In a similar fashion we can evaluate 

-3 / ^ \ ^^+1/2 



2^ '^^ 327r2z/-l 
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(209) 
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If we restrict attention to de Sitter space, then ( |207|) - (|209|) specialize to 

/$2\ ^ ^ (210) 
\ I PI 87r2 ^ ' 

+ ^ _^e-2^«(*-*°), (211) 

/ lOTT^ 

/$'2\ ^ ^ 4/fo(t-to)_ (212) 

\ Ipi 24772 V ; 

Comparison with (|166|) , (|170|) , and ( |173| ) shows that the late time results are unaf- 
fected: the dominant contribution to these expectation values comes from the infla- 
tionary sector. 

In the case of a power law inflation, we find 
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87r2' 
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^ (213) 
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^f^)f^r. (214) 




/P7 247r2 + 1 y V t / ^ ^ 

The late time constant value of ($2^ does get shifted due to (|213|) but for large p this 
shift is negligible as comparison with ( |190D shows. Even for relatively small values 
of p, this term is relatively unimportant (compare with (|192|) for p = 2). From (|198|) 
and (|204|) we know that for p ^ 1, the contributions from the inflationary sector 
to + /2 and ($'2) fall off as {to/ty and (to/^)^ respectively. These fall- 
ofTs are much slower than those given by ( |214| ) and ( pi5| ): again the pre-inflationary 
contributions are insignificant. It is only for weak power law expansions (p ~ 1) that 
this sector is of any significance. 

The rationale for the upper cutoff is simple. In quantum field theory in curved 
spacetime nontrivial ultraviolet divergences can arise because of the spacetime cur- 
vature. In principle one has to apply an appropriate regularization scheme (point- 
splitting, for example) followed by an ultraviolet subtraction. Following this more 
sophisticated procedure one obtains terms proportional to the curvature in ($^). 
However, these terms either vanish at late times (power law inflation) or are constants 
which are small compared to the contribution from the infrared sector (exponential 
inflation). In this sense our procedure is justified. (The fact that /cr] ~ 1 separates 
the low and high frequency sectors is due to the following behavior of the Hankel 
functions: oscillatory for kf] >> 1, and power law for kr] << 1.) 
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